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1.  Introduction 


Nanoparticles  have  interesting  physical  properties  that  often  vary  from  the  bulk  material.  Some 
of  these  properties,  including  increased  reactivity  {!),  are  due  to  the  high  surface  area  to  volume 
ratio  of  nanopartieles.  With  that  in  mind,  nanoparticles  may  provide  enhaneed  energy  release 
rates  for  explosive  and  propellant  reactions  (2). 

There  is  considerable  interest  in  the  self-propagating  high-temperature  synthesis  (SHS)  reactions 
of  intermetallic  compounds  because  of  the  assoeiated  energy  release  that  takes  plaee  (5)  during 
the  alloying  reaetion.  In  addition  to  the  energetic  reaction  observed  in  these  materials,  it  is 
possible  to  produce  structural  materials  that  contain  this  energy  release  property.  Onee  ignited, 
the  SHS  reaction  releases  a  large  amount  of  energy  in  a  short  period  of  time.  One  significant 
difference  between  SHS  and  typieal  eombustion  proeesses  is  that  the  reaetants  and  products  are 
confined  to  the  eondensed  state  {4).  The  SHS  proeess  has  many  potential  applications  where  heat 
generation  is  required  and  oxygen  is  not  available  or  gaseous  products  are  not  desirable.  These 
inelude  alloy  formation,  net-shape  proeessing,  propellants,  and  use  as  initiators.  One  of  the 
compounds  formed  from  the  SHS  reaction,  and  studied  here,  is  nickel  aluminide  (NiAl).  NiAl  is 
an  important  alloy  because  of  its  desirable  high  temperature  strength  and  oxidation  resistanee  (5) 
and  the  high  energy  of  formation  (6).  Reeently  Weihs  and  coworkers  have  also  used  NiAl 
nanolaminate  systems  in  applications  of  reactive  welding  (7). 

Not  surprisingly  since  the  reaction  involves  solid  starting  materials,  partiele  size  has  a  significant 
effect  on  the  properties  of  the  reaetion  product  and  the  SHS  reaction  itself  (5).  The  simulation 
and  analysis  of  nanopartiele  coalescence  without  the  SHS  reaetion  for  like  materials  is  extensive 
{9-14)  and  involves  surface  passivation  {10),  size  differences  {9,  11),  and  phase  ehange  (P) 
considerations.  The  analysis  here  ineludes  all  of  the  previously  listed  concerns  with  an  additional 
energy  release  term  from  the  heat  of  formation. 

The  foeus  of  this  report  is  to  use  atomistic  simulation  to  model  the  reactive  behavior  of  Ni-Al 
nanoparticles  in  various  eonfigurations.  Fortunately,  there  have  been  numerous  efforts  to 
determine  accurate  empirical  potentials  for  simulating  the  Ni-Al  material  system  {15).  Prior 
simulations  using  these  potentials  have  investigated  the  diffusion  of  Ni  and  A1  atoms  {15),  point- 
defect  concentrations  in  NiAl  {16),  and  plastieity  {1 7),  in  addition  to  many  other  meehanieal  and 
chemical  properties.  These  efforts  have  primarily  foeused  on  bulk  materials  rather  than 
nanoparticle  systems  {18),  even  though  there  are  many  manufacturing  processes  that  produee 
nanometer-sized  powders  for  SHS  reactions  {19).  For  this  simulation  effort  we  have  chosen  a  set 
of  embedded  atom  method  (EAM)  parameters  that  reproduce  reasonably  well  the  properties  of 
Ni,  Al,  and  NiAl  in  the  temperature  range  of  interest. 
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2.  Simulation  Approach 


In  this  work  we  employ  elassical  molecular  dynamics  (MD)  with  an  EAM  interatomic  potential 
to  study  the  SHS  reaction.  The  EAM  is  used  because  of  its  accuracy  and  capability  to  scale  up  to 
material  systems  with  over  10^  atoms.  The  MD  simulations  are  compared  with  thermodynamic 
analyses  in  order  to  provide  validation  of  the  simulation  results  and  assess  the  expected  energy 
release. 

The  MD  simulation  was  conducted  using  the  Earge-scale  Atomic/Molecular  Massively  Parallel 
Simulator  (EAMMPS)  software  package  (20).  Eor  the  Ni-Al  interactions,  the  Einnis-Sinclair 
EAM  potential  (21)  from  Angelo  et  al.  (22)  was  used.  The  Einnis-Sinclair  EAM  potential  allows 
for  non-symmetric  embedding  potential  terms,  potentially  providing  improved  accuracy  for 
metallic  alloys  (23).  In  addition  to  the  parameters  for  NiAl  from  Angelo  et  ah,  other  authors  have 
developed  parameters  for  the  Ni-Al  system  (16)  that  may  be  described  by  using  the  Einnis- 
Sinclair  EAM. 

The  three  primary  nanoparticle  sizes  that  were  considered  in  this  work  from  smallest  to  largest 
were  nanoparticles  with  1289,  5635,  and  36523  atoms  each,  which  correspond  approximately  to 
3,5,  and  10  nm,  respectively.  The  range  of  sizes  was  chosen  because  it  represents  nanoparticles 
that  may  be  produced  in  the  laboratory  and  that  offer  reasonable  computational  time  to  conduct 
parametric  studies.  Eor  the  largest  system  studied,  the  10-nm-diameter  nanoparticle  kinetic 
reaction  simulation  requires  approximately  2  days  and  64  processor  cores  to  complete  a  few 
nanoseconds  of  simulated  time  on  3.0  GHz  Intel  Woodcrest  processors. 


3.  Thermodynamic  Analysis  of  Separate  Nanoparticles 


The  separate  nanoparticle  system  is  used  as  a  model  for  powder  metallurgy  systems  where  Ni 
and  Al  particles  are  compressed  into  a  structural  component.  In  addition  to  mechanical 
properties,  the  structural  component  will  contain  stored  energy  for  future  release  through  a  SHS 
reaction.  A  thermodynamic  analysis  of  the  SHS  reaction  for  the  separate  Ni  and  Al  nanoparticle 
system  is  used  here  to  determine  the  expected  trends  and  data  points  for  simulation  validation.  In 
the  thermodynamic  analysis,  we  are  interested  in  determining  the  system  parameters  of  the  Ni-Al 
nanoparticle  system  that  contribute  to  the  combustion  temperature  and  reaction  time.  Here  we 
have  assumed  an  adiabatic  process  so  that  energy  released  to  the  surroundings  can  be  ignored. 
This  is  a  good  approximation  since  the  reaction  occurs  on  relatively  short  timescales  and  the 
nanoparticles  are  expected  to  be  included  in  a  much  larger  system  where  the  overall  surface  to 
volume  ratio  is  small,  limiting  convective  and  radiative  heat  loss.  The  validity  of  this  assumption 
is  explored  in  a  later  section.  The  SHS  reaction  of  an  equimolar  Ni  and  Al  mixture  is  written  as 
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(1) 


0.5M  +  0.5^/  ^  . 

In  order  to  compute  the  adiabatic  temperature  for  the  synthesis  reaction,  the  enthalpy  of  the 
products  and  reactants  must  be  equal: 

=  (2) 

Assuming  that  the  reaction  begins  with  the  reactants  at  600  K,  above  the  simulated  melting 
temperature  of  the  A1  nanoparticles,  the  enthalpy  of  the  reactants  is  computed  as 

Hr...  =  (0-5)(tf  )  +  (0.5)  )  =  1 1 .85  .  (3a) 

This  enthalpy  result  includes  the  enthalpy  of  solid  Ni  and  liquid  A1  (d).  The  A1  nanoparticle  is 
assumed  to  be  liquid  because  for  small  nanoparticles  the  melting  temperature  is  known  to  be 
appreciably  below  the  bulk  melting  temperature  (24).  Additionally,  for  the  EAM  potential  used 
here  (22),  the  aluminum  is  liquid  for  these  nanoparticle  sizes  at  600  K.  The  choice  of  initial 
temperature  will  have  a  nearly  linear  effect  on  the  adiabatic  temperature  as  long  as  the 
temperature  is  between  the  melting  temperature  of  the  A1  and  Ni  nanoparticles.  This  linear  effect 
has  been  observed  in  experiments  (25)  and  is  a  reasonable  assumption  so  long  as  the  heat 
capacities  of  the  solid  phases  of  Ni  and  NiAl  are  relatively  insensitive  to  temperature  in  the 
ranges  studied. 

For  the  products  of  the  SHS  process,  the  enthalpy  calculation  must  take  into  account 
contributions  from  the  melting  of  the  Ni  and  the  NiAl  nanoparticle,  enthalpy  of  formation  for  the 
NiAl  alloy,  and  changes  in  surface  energy.  The  first  of  these,  the  enthalpies  of  melting  for  Ni  and 
NiAl,  was  experimentally  determined  to  be  17.2*^„,  and  3 1.4  respectively.  The  enthalpy 

of  mixing  for  Ni  and  A1  has  garnered  close  scrutiny  in  the  experimental  community  with  a  wide 
range  of  reported  values.  The  enthalpy  of  formation  that  is  used  here  was  approximately  in  the 
middle  of  the  reported  values  at  about- 65  {6,  26,  27). 

The  last  contribution  to  the  enthalpy  of  the  products  results  from  the  change  in  surface  energy 
due  to  the  reduced  total  surface  area  of  the  combined  nanoparticle  (25).  The  contribution  to  the 
change  in  system  energy  from  the  change  in  surface  area  is  given  as  equation  4: 

^surf  ~  ^NiAI  ■  ^NiAI  ~  Al  '  ^  Al  )  '  (4) 

In  equation  4,  ,  a^,. ,  anda^,  are  the  surface  area  of  the  NiAl,  Ni,  and  A1  nanoparticles, 

respectively.  For  the  3-,  5-,  and  10-nm  A1  nanoparticles  the  reactant  surface  area  is  computed 
from  the  Gibbs  surface  {29)  as  36.32,  98.17,  and  343.7  nm  ,  respectively.  For  the  associated  Ni 
nanoparticles,  the  surface  area  is  27.15,  73.59,  and  257.87  nm^,  respectively.  The  surface  energy 
is  approximately  \\\5"y^2  for  A1  and  2573'^2  forNi  at  600  K  {20).  The  surface  area  of  the 

sintered  NiAl  nanoparticles  is  50.77,  137.18,  and  480.25  nm^  for  the  3-,  5-,  and  10-nm 
nanoparticle  cases,  respectively.  In  experimental  analysis  of  the  free  surface  energy  of  NiAl  near 
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its  melting  point,  the  free  surface  energy  has  been  reported  as  1400  ""^2  {SI).  The  approximate 

change  in  energy  versus  nanoparticle  size  is  tabulated  in  table  1.  In  table  1,  the  trend  is  for  a 
lower  surface  energy  contribution  to  the  reaction  as  the  nanoparticle  size  increases.  Intuitively, 
one  may  expect  this  because  the  surface  area  to  volume  ratio  is  also  decreasing  with  increasing 
particle  size,  and  therefore  has  less  influence  on  the  sintering  process.  With  the  enthalpy  of 
formation  for  NiAl  around  -  65*%,^, ,  the  surface  energy  contribution  to  the  change  in  enthalpy 

for  coalescence  of  10-nm-diameter  nanoparticles  is  less  than  10%  of  the  total  enthalpy  change. 
This  means  that,  even  at  relatively  small  nanoparticle  sizes,  e.g.,  10  nm,  the  effect  of 
nanoparticle  size  on  energy  release  is  minimal. 


Table  1.  Change  in  surfaee  energy  versus  nanopartiele  size. 


Nanopartiele  Radius 
(nm) 

(kJ/mol) 

3 

-18.35 

5 

-11.41 

10 

-6.17 

With  the  preceding  discussion,  it  is  possible  to  take  into  account  many  of  the  sources  of  enthalpy 
change  in  the  reaction  products  including  phase  and  surface  area  changes.  The  enthalpy  of  the 
products  is  now  estimated  as 

'Taa 

^  prod  ~  ^  form, NiAl  surf  ^  ^  ^  melt,Ni  '  (6) 

29SK 

The  heat  capacity  for  solid  and  liquid  NiAl  is  given  in  Kubaschewski  et  al.  (21).  For  the  3-nm 
case,  assuming  the  NiAl  nanopartiele  melting  temperature  to  be  about  1350  K,  or  the  melting 
point  of  a  similarly  sized  Ni  nanopartiele,  it  is  possible  to  compute  the  adiabatic  reaction 
temperature  (table  2).  Notice  in  table  2  that  if  no  surface  energy  contribution  is  considered,  i.e., 
infinitely  large  spheres,  the  final  adiabatic  temperature  is  computed  to  be  1599  K.  In  section  6, 
we  will  observe  that  these  results  are  reasonable  and  accurately  predict  the  simulated  increase  in 
temperature  attributable  to  the  contribution  from  the  surface  energy. 

Table  2.  Computed  adiabatic  temperature  versus  uauoparticle  radius, 
iucludiug  coutact  of  flat  surfaces  or  iufiuitely  sized  spheres. 


Nanopartiele  Radius 

Tad 

(nm) 

(K) 

3 

2115 

5 

1920 

10 

1772 

00 

1599 

4 


4.  The  Coalescence  Processes 


For  Ni  and  A1  nanoparticles,  the  SHS  reaction  consists  of  two  processes,  namely,  eoalescence 
and  alloying.  In  this  work,  we  have  considered  the  eoalescence  of  a  two  nanoparticle  system  with 
an  A1  and  a  Ni  nanopartiele  with  an  atomic  ratio  of  unity.  A  complete  SFIS  reaction  of  this 
system  will  result  in  a  single  NiAl  nanoparticle.  The  MD  simulations  used  to  work  model 
adiabatic  conditions  with  constant  number  of  atoms  and  total  system  energy.  The  purpose  of 
these  simulations  is  to  analyze  the  effeet  of  nanopartiele  size  on  sintering  time  and  adiabatic 
combustion  temperature,  and  to  visualize  the  proeess.  The  assumed  proeess  is  illustrated  in 
figure  1 .  In  figure  I ,  the  nanopartieles  are  initially  in  contact  at  a  point  (a)  and  the  A1 
nanopartiele  is  larger  than  the  Ni  nanopartiele  because  of  the  longer  Al-Al  bond  length.  The 
simulations  are  initialized  at  600  K  so  that  the  A1  nanopartiele  is  liquid  and  the  Ni  nanopartiele  is 
solid.  In  figure  I,  the  sintering  process  proeeeds  with  the  liquid  A1  nanopartiele  initially  coating 
the  solid  Ni  nanopartiele  while  forming  some  Ni-Al  bonds  on  the  surfaee  (b-d).  Next,  the 
alloying  proeess  proeeeds  with  the  Ni  nanopartiele  being  heated  above  its  melting  point  and 
becoming  liquid  so  that  mixing  may  occur  (e).  The  formation  of  Ni-Al  bonds  beyond  the 
interfacial  surface  requires  diffusion  of  A1  into  the  Ni  nanopartiele  or  Ni  into  the  liquid  Al.  Either 
of  these  proeesses  is  possible  but  since  diffusion  is  a  relatively  slow  proeess  in  solid  materials  it 
is  expected  that  the  Ni  nanopartiele  must  melt  before  the  eoalescenee  proeess  proeeeds 
appreeiably. 


Figure  1.  Illustration  of  sintering  proeess  showing  liquid  Al  nanopartiele  first  eoating  the  solid  Ni  nanopartiele  and 
then  eomplete  alloying  after  the  Ni  nanopartiele  has  melted. 

The  nanopartiele  sintering  proeess  is  driven  by  two  sources  of  energy  as  previously  diseussed. 
The  first  of  these  is  a  deerease  in  surface  area  that  lowers  the  total  surface  energy  of  the  system. 
This  energy  release  meehanism  is  also  observed  in  the  sintering  of  homogeneous  material 
systems  sueh  as  silicon  nanoparticles  {28,  32).  The  second  source  of  energy  is  from  the  reactive 
synthesis  that  oeeurs  initially  at  the  interfaee  between  the  nanoparticles  and  later  throughout  the 
entire  system.  The  energy  release  from  the  surfaee  sintering  is  proportional  to  the  surface  area  of 
the  Ni  nanopartiele  that  is  eoated  by  Al  and  in  the  whole  system  to  the  total  number  of  Ni  and  Al 
atoms.  Additionally,  with  the  temperature  increase  there  is  a  decrease  in  the  viseosity  of  the 
liquid  Al  that  will  affeet  the  predicted  coalescence  time. 
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The  coalescence  of  nanoparticles  in  the  liquid  and  solid  phases  has  been  examined  extensively 
{9-11).  These  studies  are  primarily  concerned  with  the  coalescence  of  two  liquid  or  two  solid 
nanoparticles.  The  analysis  for  the  Ni-Al  system  requires  considering  the  coalescence  of  a  liquid 
A1  nanoparticle  and  a  solid  Ni  nanoparticle.  Lewis  et  al.  (P)  considered  the  coalescence  of  a 
liquid  and  a  solid  gold  nanoparticle;  this  is  similar  to  the  situation  here  except  that  the  material 
system  considered  was  homogeneous. 

In  Lewis  et  al.  (P),  the  author  was  able  to  simulate  two  phases  occurring  simultaneously  for  a 
single  material  by  choosing  the  size  of  each  nanoparticle  such  that  at  a  specific  temperature  the 
phase  of  the  nanoparticles  was  different.  Lewis  found  that  coalescence  proceeded  in  two  stages: 
first,  the  contact  area  was  maximized;  and  second,  “sphericization”  took  place  driven  by  surface 
diffusion.  The  first  stage  is  much  faster  than  the  second  and  is  very  similar  to  the  process 
observed  here,  where  the  Al  nanoparticle  maximizes  the  contact  area  and  partially  coats  the  Ni 
nanoparticle.  In  this  case,  there  is  an  added  driving  force  other  than  the  surface  energy, 
specifically  the  energy  release  on  forming  of  Ni-Al  bonds  as  compared  to  the  Al-Al  and  Ni-Ni 
bonds.  During  the  second  stage,  the  atoms  in  the  two  nanoparticles  diffuse  and  rearrange  until 
the  system  becomes  a  single  spherical  nanoparticle.  This  stage  is  driven  strongly  by  the 
formation  of  Ni-Al  bonds  and  is  expected  to  occur  on  a  much  shorter  timescale  than  for  two 
nanoparticles  of  the  same  material.  The  analytical  model  and  MD  simulation  results  shown  in 
sections  5  and  6  explore  this  assumption. 


5.  Phenomenological  Model  of  Nanoparticle  Reactive  Sintering 


To  gain  further  insight,  we  have  developed  a  phenomenological  model  for  the  reactive  sintering 
of  Ni  and  Al  nanoparticles.  The  model  includes  energy  release  from  surface  energy,  bond 
formation,  and  viscous  dissipation  through  deformation.  Frenkel  (33)  has  developed  a  model  for 
the  coalescence  of  two  homogeneous  nanoparticles;  however,  his  model  did  not  account  for  any 
phase  change,  kinetic  sintering,  or  heterogeneous  materials.  Here  we  extended  Frenkel’s  model 
to  consider  the  coalescence  of  two  liquid  drops  to  take  into  account  the  coalescence  of  a  liquid 
and  a  solid  drop  with  reactive  synthesis. 

The  analytical  model  is  initialized  with  the  Al  and  Ni  nanoparticles  in  contact  at  a  point.  The 
distance  from  nanoparticle  center  to  center  is  equal  to  the  sum  of  the  respective  radii,  denoted  as 
D  in  figure  2.  The  sintering  process  initially  proceeds  by  the  liquid  Al  nanoparticle  coating  the 
solid  Ni  nanoparticle,  as  illustrated  in  figure  1.  During  this  phase  of  the  sintering  process,  two 
sources  of  energy  release  are  occurring;  the  first  of  these  is  related  to  the  decrease  in  surface  area 
and  is  proportional  to  the  respective  surface  tension  values.  The  second  source  of  energy  release 
is  from  the  formation  of  Ni-Al  bonds  at  the  interfacial  region.  Figure  3  is  an  illustration  of  the 
parameters  used  to  model  the  coalescence  time. 
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Figure  2.  Illustration  of  parameters  used  in 

analytical  model  of  reactive  coalescence 
of  Ni  and  A1  nanoparticles. 


(b) 


Figure  3.  (a)  Plot  of  exposed  Ni  and  A1  nanoparticle  surface  area  as  a  function  of  distance  between  nanoparticle 

centers  and  (b)  a  plot  of  total  exposed  surface  area  as  a  function  of  distance  between  nanoparticle  centers. 
These  results  assume  a  Ni  nanoparticle  of  radius  4.53  nm  and  an  A1  nanoparticle  of  5.23  nm.  Notice  that 
the  total  exposed  surface  area  is  monotonically  decreasing,  indicating  that  the  surface  energy  is  also 
decreasing  monotonically 


In  figure  2,  2a  is  the  diameter  of  a  circle  circumscribed  by  the  contact  circumference  of  the  two 
nanoparticles.  Vai  and  vni  are  the  distance  from  the  A1  and  Ni  nanoparticle  surface  to  the  surface 
of  the  contact  circle,  respectively.  0  is  the  contact  angle  as  measured  from  the  center  of  the  Ni 
nanoparticle  and  ranges  from  0  to  ti  radians.  In  order  to  model  the  change  in  energy  of  the 
coalescing  nanoparticle  system,  three  energy  change  mechanisms  must  be  considered.  These 
mechanisms  are  energy  release  due  to  change  in  surface  area,  energy  release  due  to  kinetic 
reactions  at  the  interface,  and  energy  loss  due  to  viscous  dissipation.  The  rate  of  energy  change 
due  to  all  three  must  balance  at  all  times. 

The  first  energy  term  considered,  namely,  the  surface  energy  of  the  nanoparticle  system,  is 
simply  the  surface  tension  times  the  total  exposed  surface  area.  This  energy  term  is  written  as  a 
sum  of  the  A1  and  Ni  nanoparticle  contributions: 
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^ surf  ^ Ni  ^ Ni^ex^iosed  ^ Al  ^ Al , exposed  * 

(10) 

The  exposed  area  of  the  Ni  nanopartiele  can  be  written  as 

^Ni,exposed  ~  ^Ni  ~  ~  ~  ’ 

(11) 

where 

V*  =^m(1-cos(0)). 

(12) 

Initially,  during  the  sintering  process,  the  Ni  nanopartiele  is  assumed  to  remain  in  the  solid 
phase,  thus  maintaining  a  constant  radius.  This  assumption  is  reasonable  because  of  the  higher 
melting  temperature  of  the  Ni  nanopartiele. 

The  exposed  surface  area  of  the  A1  nanopartiele  is  written  as 

Sai, exposed  =  ^.4/  “  -  2nr^,v , 

(13) 

where 

CN 

1 

1 

II 

(14a) 

(14b) 

The  radius  of  the  A1  nanopartiele  is  eomputed  numerieally  by  using  eonservation  of  volume  for 
the  A1  nanopartiele.  The  exposed  surfaee  area  of  eaeh  nanopartiele  versus  the  eenter  to  eenter 
distanee  is  plotted  in  figure  3a.  Notiee  that  although  the  surfaee  area  of  the  A1  nanopartiele 
inereases  during  most  of  the  eoaleseenee  proeess  the  eombined  total  surfaee  area  of  the  Ni  and 
A1  nanopartieles  deereases  monotonieally  throughout  the  entire  eoaleseenee  proeess.  In  figure  3, 
the  eenter-to-eenter  distanee  never  reaehes  zero  beeause  the  eoaleseenee  is  eonsidered  eomplete 
onee  the  Ni  nanopartiele  is  eompletely  enveloped  by  the  A1  nanopartiele. 

The  seeond  souree  of  energy  release,  namely,  the  reaetive  synthesis  term,  is  eonsidered  by 
assuming  a  eonstant  surfaee  density  of  the  Ni  nanopartiele  and  the  transient  eontaet  area  of  the 
Ni-Al  interfaee: 


^reactive  ~  Pm  ,surface  ^  int  erface 


K 


hondenergy  * 


(15) 


The  surfaee  density  term,  ,  is  proportional  to  the  number  of  Ni-Al  bonds  at  the  eontaet 

interfaee.  The  surfaee  density  and  bond  energy  terms,  ,  can  be  combined  into  a  single 

constant  that  defines  the  energy  release  per  unit  area  of  interface: 


density  P Ni,surface  ^bondenergy  * 


(16) 
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The  interfacial  contact  area  is  a  function  of  the  distance  between  nanoparticle  centers  (figure  4). 
The  interfacial  area  increases  monotonically  until  the  Ni  nanoparticle  surface  is  completely 
covered.  This  result  is  expected  since  the  reactive  energy  term  is  negative,  or  releases  energy 
during  the  entire  process,  in  addition  to  the  minimization  of  surface  energy  that  is  driven  by  the 
surface  tension  of  Ni  and  Al.  The  interfacial  area  is  written  as 


^mterface  ~  Ni^Ni  ’ 

where  is  a  function  of  ^as  given  in  equation  12. 


(17) 


Figure  4.  Contact  or  interface  area  as  a  function  of  center-to-center 
distance.  The  contact  area  is  increasing  as  the  nanoparticles 
move  closer  together  (right  to  left  on  x-axis). 


The  third  energy  term  represents  the  viscous  dissipation  due  to  deformation  of  the  Al 
nanoparticle.  This  viscous  dissipation  is  a  function  of  the  viscosity  in  the  liquid  Al  nanoparticle 
and  the  rate  of  deformation.  The  extent  of  the  viscous  flow  can  be  specified  by  the  decrease  in 
distance  between  the  center  of  each  drop,  and  the  surface  of  contact  with  the  Ni  nanoparticle.  A 
velocity  gradient,  / ,  can  be  defined  as(^T))/r^, .  The  energy  dissipated  in  the  whole  body  per  unit 

time  is,  therefore,  approximately 


dE„ 


dt 


^AI,Q  / 

2r]  =  ■ 

0  V 


dt  y 


(18) 


where  r]  is  the  viscosity  of  liquid  Al  and  ^  is  the  initial  radius  of  the  Al  nanoparticle. 


By  conservation  of  energy,  the  rate  of  coalescence  can  now  be  computed: 
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dE... 


dt 


dt  dt 


s 

3 


'Al,od\ 


dD 

dt 


\2 


d_ 

dt 


■S 


M,exp  osed 


+  (7 


I  .  e 

Al  Al, exposed 


V-\m 

‘  dr 


(19a) 

(19b) 


After  writing  equation  19b  in  terms  of  d0/dt  and  simplifying  the  right  and  left  hand  sides,  we 
find  that  equation  19b  is  only  linearly  dependent  on  d0/dt.  Even  with  this  simplification, 
equation  19  is  most  easily  solved  numerically  using  an  iterative  solver.  In  order  to  solve  equation 
19b,  we  need  some  physical  properties  of  Al,  Ni,  and  NiAl.  The  dynamic  viscosity  of  bulk 
molten  Al  at  the  melting  temperature  is  about  rj  =  1.3- 10“^ Pa -5  (54).  Based  upon  a  comparison  of 

the  configurational  energy  in  MD  simulations  of  separate  nanoparticles  and  Al-coated  Ni 
nanoparticles,  the  energy  release  per  unit  area,  ,  is  estimated  to  be  20.7«(i(^2 .  This  number 

is  computed  by  subtracting  the  system  energy  of  an  Al-coated  Ni  nanoparticle  system  from  the 
energy  of  a  system  with  separate  nanoparticles  and  dividing  by  the  interfacial  surface  area.  This 
method  results  in  the  net  change  in  energy  during  coating  of  the  Ni  surface  with  Al  since  some 
Al-Al  bonds  are  lost  during  the  coating  process  while  some  Ni-Al  bonds  are  formed  at  the 
interface.  By  numerically  solving  equation  19b,  we  are  able  to  compute  the  contact  angle,  0,  as  a 
function  of  time  and  relate  this  to  total  exposed  surface  area  of  the  coalescing  nanoparticles.  This 
result  is  presented  in  figure  5  along  with  a  comparison  to  the  MD  simulation  results. 


Figure  5.  Total  system  surface  area  versus  time  from  mathematical 

model  and  MD  simulations  for  the  sintering  of  10-nm-diameter 
nanoparticles,  where  the  final  surface  area  of  the  NiAl 
nanoparticle  is  approximately  480  nm^. 
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Although  qualitatively  the  results  in  figure  5  show  similar  trends,  the  absolute  rate  of 
coalescence  is  slightly  under  predicted  by  the  model.  This  difference  can  be  attributed  to  the 
obvious  simplicity  of  the  model  and  more  specifically  to  the  difficulty  in  obtaining  accurate 
material  parameters.  For  instance,  it  is  difficult  to  compare  the  viscosity  of  a  nanoparticle  to  the 
bulk  material  {12),  and  since  the  coalescence  time  is  linearly  dependent  upon  the  viscosity  a 
change  in  viscosity  is  directly  proportional  to  a  change  in  modeled  coalescence  time. 
Additionally,  the  energy  release  per  unit  area  term  assumes  that  the  net  change  in  energy  due  to 
the  addition  of  Ni-Al  bonds  at  the  interface  is  a  constant  value.  This  is  likely  not  completely 
accurate  since  fewer  A1  bonds  must  be  broken  to  form  new  Ni-Al  bonds  during  the  initial  contact 
of  the  nanoparticles.  Fiowever,  the  deviation  in  this  energy  release  term  is  likely  to  be  minimal. 
The  deviation  of  the  model  time  from  the  simulation  results  at  about  50  ps  is  due  to  the  switch 
from  stage  1  to  stage  2  in  the  kinetic  coalescence  process  as  described  by  Lewis  et  al.  (P).  As 
described  by  Lewis,  during  stage  2,  surface  diffusion  is  the  predominant  factor  in  continued 
coalescence  and  is  a  much  slower  process  than  contact  area  maximization.  The  actual  simulation 
results,  as  compared  with  the  illustration  in  figure  1,  of  the  observed  coalescence  process  are 
given  in  figure  6. 


Initial  Stage  1  Stage  2  Completed 

Fast  (50  ps)  Slow  (450  ps)  Slow  (400+  ps) 


Figure  6.  Cross-sectional  view  from  MD  simulations  of  Ni/Al  nanoparticle  sintering  process  showing  the  start  of 
the  second  stage  of  coalescence  where  diffusion  is  the  driving  force  as  opposed  to  contact  area 
maximization.  Al  atoms  are  blue  and  Ni  atoms  are  red. 

In  figure  6,  each  of  the  steps  in  the  coalescence  process  is  shown  with  plots  from  an  MD 
simulation  of  the  coalescence  of  10-nm-diameter  Al  and  Ni  nanoparticles.  The  correlation  of  the 
sintering  stages  to  the  reaction  temperature  and  time  is  illustrated  in  figure  7  for  the  sintering  of 
separate  10-nm-diameter  nanoparticles.  In  the  initial  step,  the  liquid  Al  nanoparticle  (the  blue 
atoms  in  figure  6)  has  melted  and  is  spherical  in  shape.  The  solid  Ni  nanoparticle  (the  red  atoms) 
has  large  faceted  sides  and  is  a  single  crystal,  a  typical  configuration  for  a  crystalline 
nanoparticle  at  low  temperatures.  During  stage  1,  the  Al  nanoparticle  is  attracted  to  the  Ni 
surface  because  of  the  dual  driving  forces  of  surface  energy  minimization  and  Ni-Al  bond 
formation.  This  period  lasts  about  50  ps  in  this  simulation  as  noted  in  figures  6  and  7.  Between 
stages  1  and  2,  the  driving  forces  associated  with  the  surface  energy  are  counteracted  by  a 
resistance  to  flow  in  the  Al  nanoparticle,  causing  the  coalescence  process  to  slow  down 
dramatically.  During  stage  2,  lasting  about  450  ps,  the  surface  area  is  not  changing  so  that  energy 
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release  from  the  surface  energy  terms  has  ceased  to  contribute  to  the  change  in  system  potential 
energy.  The  subsequent  energy  release  is  entirely  attributable  to  the  formation  of  Ni-Al  bonds. 
This  stage  lasts  a  much  longer  time  than  the  initial  nanoparticle  coalescence  stage  and  is 
governed  by  the  material  diffusion  coefficients.  Initially,  at  stage  2,  the  Ni  nanoparticle  is  still 
solid  and  the  formation  of  Ni-Al  bonds  is  only  possible  by  A1  diffusing  into  the  Ni  core  or  Ni  on 
the  surface  of  the  core  melting  and  diffusing  away  from  the  interface.  This  process  proceeds  until 
the  Ni  core  has  reached  its  melting  point  and  mixing  of  the  remaining  Ni  and  A1  atoms  occurs 
more  rapidly,  driven  by  the  enthalpy  of  formation  of  NiAl.  From  stage  2  until  complete  alloying 
has  occurred,  taking  approximately  400  ps,  diffusion  and  mixing  of  Ni  and  A1  atoms  is  the 
primary  driving  force. 


Figure  7.  Time  vs.  temperature  plot  for  sintering  of  separate 

10-nm-diameter  A1  and  Ni  nanoparticles.  The  various 
stages  if  the  coalescence  processes  are  denoted  on  the 
curve,  including  the  final  completion  stage  that  occurs 
after  the  Ni  nanoparticle  has  melted. 


6.  MD  Simulation  Results  of  Separate  Nanoparticle  Reactivity 


We  have  previously  predicted  the  adiabatic  temperature  and  sintering  time  for  the  reactive 
sintering  process  of  separate  equimolar  nanoparticles  of  A1  and  Ni.  In  figure  8,  the  MD 
simulation  results  for  the  equimolar  nanoparticles  are  plotted  along  with  the  computed  adiabatic 
temperature  for  each  considered  particle  size. 
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Figure  8.  Temperature  vs.  time  in  the  sintering  of  nanoparticles 
with  an  Ni:Al  ratio  of  1 : 1 . 

Note:  The  subscripts  in  the  legend  refer  to  the  number  of  atoms  of  each 
material  and  correspond  to  nanoparticles  of  diameter 
approximately  3,  5,  and  10  nm.  The  color-coded  dashed  lines  are 
the  computed  adiabatic  temperature  from  the  thermodynamic 
analysis.  The  black  dashed  line  is  the  predicted  temperature  for 
coalescence  of  bulk  A1  and  Ni. 

From  figure  8,  it  is  apparent  that  the  predicted  adiabatic  temperature  is  in  close  agreement  with 
the  simulated  temperature.  Variability  of  the  computed  temperature  arises  from  the  wide  range  of 
experimental  results  for  the  surface  tension  for  liquid  A1  and  solid  Ni,  the  reported  enthalpy  of 
formation  for  NiAl,  and  the  assumed  melting  temperature  for  the  Ni  and  NiAl  materials  at  this 
scale.  Each  of  these  experimental  data  points  is  used  in  the  thermodynamic  analysis  and 
contributes  to  the  small  inaccuracies  in  the  predicted  temperature. 

The  characteristic  time  for  reactive  synthesis  that  we  use  here  is  defined  by  Zhao  et  al.  (35)  as  t 
when 


r(t)  =  ?;+ 0.8(7; -rj,  (20) 

where  is  the  initial  temperature,  is  the  maximum  size  dependent  temperature  reached,  and 
r(t)is  the  transient  temperature.  The  computed  reaction  times  are  given  in  figure  9  and  illustrate 
that  the  time  required  for  separate  nanoparticles  to  react  has  a  power  law  relationship  between 
nanoparticle  volume  (3)  and  surface  area  (2).  This  implies  that  not  only  will  the  reaction 
temperature  be  higher,  but  it  will  occur  more  rapidly  with  decreases  in  particle  size,  to  a  power 
of  about  2.5.  This  is  important  because  a  high  rate  of  energy  release  is  desirable  for  many 
applications. 
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Figure  9.  Reaction  time  versus  (a)  A1  nanoparticle  diameter,  (b)  A1  nanoparticle  surface  area,  and  (c)  number  of  A1 

atoms.  Note  the  nearly  linear  relationship  (dashed  line)  of  reaction  time  between  number  of  atoms  (volume) 
and  surface  area. 


By  observing  the  MD  simulations  and  analyzing  the  shape  of  the  eurves  in  figure  8  for 
temperature  versus  time,  we  have  surmised  that  there  are  two  reaetion  rates  to  eonsider.  The  first 
is  during  the  eoalescenee  proeess  (ealled  the  growth  rate  [see  0  to  500  ps  for  AI36523  curve  in 
figure  8])  and  Ni  nanoparticle  melting,  and  the  second  is  the  rapid  formation  of  the  NiAl  alloy 
from  liquid  A1  and  Ni  (convergence  rate  [see  500  ps  and  later  for  the  AI36523  curve  in  figure  8]). 
The  temperature  at  which  the  transition  occurs  is  size  dependent  because  the  melting  temperature 
of  the  Ni  nanoparticle  is  also  size  dependent.  If  the  process  is  not  perfectly  adiabatic,  and  some 
heat  is  lost  to  the  surroundings,  it  is  also  possible  that  the  first  process  would  not  proceed  far 
enough  for  the  Ni  nanoparticle  to  melt,  and  thus  the  reaction  would  halt.  This  would  only  occur 
with  larger  nanoparticles  that  require  longer  reaction  times  during  which  some  energy  loss  to  the 
surroundings  is  likely.  This  is  an  important  consideration  in  real-world  applications  that  are  not 
perfectly  adiabatic,  but  when  complete  alloying  is  desired.  One  reason  for  this  observed  increase 
in  reaction  rate  is  because  the  heat  generated  from  the  formation  of  Ni-Al  bonds  will  conduct 
into  the  core  nanoparticle  so  that  when  the  reaction  front  reaches  the  inner  atoms  they  will  have  a 
higher  diffusion  coefficient,  which  in  turn  increases  the  reaction  rate. 


7.  Reactive  Sintering  of  Core-Shell  Nanoparticles 


7.1  Aluminum-coated  Nickel 

In  this  section,  we  discuss  the  sintering  process  for  an  Al-coated  Ni  nanoparticle  followed  by  a 
discussion  of  a  Ni-coated  A1  nanoparticle.  Both  of  these  systems  can  be  used  as  a  model  for 
highly  compacted  Ni  and  A1  nanoparticles  or  one  material  serving  as  a  matrix  for  nanoparticles 
of  the  other.  In  the  first  model  system,  we  assume  that  a  Ni  nanoparticle  has  been  coated  with  A1 
and  equilibrated  without  the  Ni  melting  or  any  further  reaction  occurring.  Results  for  the  reaction 
time  and  temperature  are  presented  and  a  comparison  with  the  separate  nanoparticle  case  is 
given.  Here  again,  we  have  considered  three  system  sizes  with  1289,  5635,  and  36523  atoms 
each  of  A1  and  Ni. 
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An  initial  estimate  is  that  the  coalescence  process  for  the  fully  coated  nanoparticle  system  will  be 
a  truncated  version  of  the  separate  nanoparticle  case.  In  the  coated  nanoparticle  system,  we  do 
not  have  the  first  stage  of  coalescence  occurring  and  only  observe  the  second  stage,  namely, 
diffusion  of  Ni  and  A1  atoms  to  form  Ni-Al  bonds.  The  sintering  temperature  versus  time  plot  is 
given  in  figure  10  and  shows  an  interesting  result.  Whereas  the  maximum  temperature  reached 
increases  with  decreasing  nanoparticle  size  for  coalescence  of  separate  nanoparticles,  the 
opposite  is  true  here;  the  temperature  decreases  with  decreasing  nanoparticle  size. 


Figure  10.  Temperature  vs.  time  in  the  sintering  of  Al-eoated 
Ni  nanopartieles  with  an  Ni:Al  ratio  of  1 : 1 . 


In  figure  10,  the  observed  decrease  in  adiabatic  temperature  is  due  to  the  fact  that  the  ratio  of 
atoms  near  the  interfacial  region  to  the  atoms  in  the  bulk  nanoparticle  decreases  as  the 
nanoparticle  size  increases.  Atoms  in  the  interfacial  region  have  already  formed  Ni-Al  bonds  and 
are,  therefore,  already  at  a  lower  configurational  energy  than  if  they  were  contained  in  a 
homogeneous  nanoparticle  of  either  pure  A1  or  Ni.  If  we  extend  the  adiabatic  temperature 
relationship  to  infinitely  large  particles,  we  would  approach  the  result  obtained  from  the  analysis 
of  separate  nanoparticle  as  they  increase  in  size.  The  thermodynamic  analysis  is  similar  except 
that  the  surface  energy  term  is  zero  and  the  enthalpy  of  formation  is  lowered  by  a  factor 
proportional  to  the  ratio  of  surface  area  to  volume.  The  enthalpy  of  the  products,  equation  6, 
modified  for  coated  nanopartieles  becomes 
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where  t  is  a  computed  thickness  value  for  the  interfacial  layer,  is  the  area  of  the  interfacial 

region,  and  V  is  the  volume  of  the  Ni  core.  In  order  to  determine  the  correct  empirical  thickness 
value,  t,  for  equation  21  we  have  used  the  adiabatic  temperature  computed  in  the  MD  simulation 
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results  for  the  Al-eoated  Ni  nanopartiele.  These  results  indieate  that  an  interfaee  thiekness  of 
0.07  nm  is  able  to  aeeurately  prediet  the  adiabatie  temperature  observed  in  the  MD  simulations 
(figure  1 1).  In  figure  1 1,  it  is  apparent  that  the  adiabatie  eombustion  temperature  is  highly  size 
dependent  for  nanopartieles  of  less  than  10  nm  in  diameter.  For  very  small  nanopartieles,  less 
than  1  nm  diameter,  there  is  little  predieted  ehange  in  temperature  from  the  initial  temperature  of 
600  K  sinee  most  of  the  potential  Ni-Al  bonds  have  already  been  formed. 
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Figure  1 1 .  Comparison  of  thermodynamically  determined 

adiabatic  temperature  for  Al-coated  Ni  nanopartiele 
and  results  from  MD  simulation. 

From  the  sintering  of  separate  nanopartieles,  it  is  expeeted  that  the  reaetion  time  will  be  linearly 
related  to  the  radius  of  the  nanopartiele  to  a  power  of  about  2.5.  In  figure  12,  this  appears  to  be 
the  ease  for  this  range  of  nanopartiele  sizes.  A  slight  deviation  from  the  separate  nanopartiele 
result  is  probably  related  to  the  faet  that  the  eoaleseenee  proeess,  stage  1,  is  not  ineluded  in  this 
model  system  and  diffusion  takes  longer  to  initialize  the  kinetie  reaetion  proeess. 
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Figure  12.  Reaction  time  versus  number  of  A1  atoms  in  the 
Al-coated  Ni  nanoparticle  system. 

The  results  for  the  Al-eoated  Ni  nanopartiele  indieate  the  trends  that  one  might  expeet  from  a 
material  system  that  ineluded  an  A1  matrix  with  embedded  Ni  nanopartieles.  From  the  results  in 
figures  1 1  and  12,  there  are  two  eompeting  reaetion  results,  namely,  reaetion  time  and  maximum 
temperature.  In  figure  12,  we  see  that  as  the  Ni  nanopartiele  size  deereases  the  reaetion  time 
deereases,  eausing  the  energy  release  rate  to  inerease.  A  seeond  observation  that  ean  be  made 
from  figure  1 1  is  that  the  reaetion  temperature  deereases  with  deereasing  Ni  nanopartiele  size, 
potentially  minimizing  the  effeet  of  the  rapid  energy  release. 

Looking  more  elosely  at  the  reaetion  time  versus  number  of  atoms  for  the  separate  nanopartiele 
and  Al-eoated  Ni  nanopartiele  eases,  we  observe  a  similar  relationship  of  reaetion  time  to 
nanopartiele  size  as  that  found  in  separate  nanopartieles.  In  both  oases,  the  reaetion  time  appears 
to  have  a  power  law  relationship  with  radius,  with  an  exponent  of  2.5.  The  aooelerated 
temperature  inerease  in  figure  10  after  about  900  ps  for  the  AI36523  curve  is  the  oonvergenoe  rate 
discussed  previously. 

7.2  Nickel-coated  Aluminum 

The  Ni-ooated  A1  nanopartiele  system  has  garnered  some  interest  beoause  when  A1  melts  there  is 
an  experimentally  observed  inerease  in  volume  of  about  6.5%.  This  inerease  in  volume  oreates  a 
large  stress  in  the  Ni  ooating  and  may  result  in  oatastrophio  failure  and  fragmentation  of  the 
nanopartiele  (18).  In  the  work  by  Delogu  (18),  the  fragmentation  is  only  observed  for  a  speeifie 
set  of  eonditions  ineluding  Ni  shell  thiekness.  In  this  work,  even  with  the  appropriate  Ni  shell 
thiekness,  fragmentation  of  the  nanopartiele  is  not  observed.  With  the  empirieal  potentials  and 
parameters  eonsidered  here,  the  kinetie  alloying  reaetion  is  fast  and  results  in  a  liquid  NiAl 
nanopartiele  without  fragmentation. 


17 


In  the  first  set  of  analyses,  the  Ni  shell  contains  the  same  number  of  atoms  as  the  A1  core.  This 
results  in  an  A1  core  with  a  radius  of  2.8  nm  and  a  Ni  shell  of  about  0.7  nm.  Since  the  bond 
length  of  Ni  is  less  than  Al,  the  shell  is  thinner  and  the  contact  area  is  initially  greater  than 
observed  in  similarly  sized  Al-coated  Ni  nanoparticles.  As  with  the  nanoparticle  coalescence 
simulations,  the  temperature  of  the  system  is  raised  to  600  K,  at  which  time  a  constant  energy 
simulation  is  used  to  analyze  the  energy  conversion  rate  and  the  adiabatic  temperature  rise  of  the 
system. 

In  the  initial  simulations  with  an  atomic  ratio  of  unity,  there  are  more  Ni-Al  bonds  in  the  Ni- 
coated  nanoparticle  than  the  Al-coated  nanoparticle.  It  may,  therefore,  be  expected  that  the  total 
system  energy  would  initially  be  lower  in  the  Ni-coated  Al  nanoparticle  system  than  the  Al- 
coated  Ni  system.  In  fact,  the  opposite  is  true,  because  although  there  are  more  Ni-Al  bonds  in 
the  Ni-coated  system  there  are  fewer  Ni-Ni  bonds  than  in  the  Al-coated  nanoparticle.  Since  Ni- 
Ni  bonds  are  stronger  than  Al-Al  bonds,  the  total  initial  energy  is  lower  in  the  Al-coated  Ni 
nanoparticle  system.  Since  the  final  configuration  of  both  systems  is  a  completely  alloyed  NiAl 
nanoparticle,  the  system  energy  change  for  the  Ni-coated  nanoparticle  is  greater  than  the  Al- 
coated  system.  This  greater  change  in  potential  energy  results  in  the  computed  adiabatic 
temperature  for  the  Ni-coated  Al  nanoparticle  being  higher  than  the  Al-coated  nanoparticle 
(figure  13),  although  the  difference  is  not  large. 


Figure  13.  Comparison  of  combustion  temperature  for  Ni-coated 
Al  nanoparticle  versus  Al-coated  Ni  nanoparticle. 

In  the  two  following  simulations,  we  simulated  an  approximately  5-nm-diameter  Al  nanoparticle 
coated  with  either  a  1-  or  2-nm-thick  Ni  coating.  For  these  simulations  the  temperature  was 
controlled  using  an  NVT  ensemble.  This  temperature  control  was  used  in  order  to  rapidly 
increase  the  temperature  from  about  300  K  to  above  the  melting  point  of  the  Al  nanoparticle  at  a 
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rate  of  0.1  K/ps.  When  the  A1  nanopartiele  melts  and  expands,  a  large  sudden  inerease  of  stress 
in  the  Ni  eoating  is  observed.  For  the  2-nm-thiek  Ni  eoating,  the  inerease  in  stress  is  not  high 
enough  to  cause  failure  of  the  coating.  In  this  case,  the  outer  shell  expands  slightly  but  does  not 
crack.  For  the  1-nm-thick  case,  the  stress  in  the  Ni  shell  is  high  enough  to  cause  failure.  When 
the  shell  fails,  the  A1  begins  to  leak  out  onto  the  surface  of  the  nanopartiele  but  no  fragmentation 
is  observed.  If  the  nanopartiele  were  surrounded  by  oxygen,  this  may  initiate  the  A1  oxidation 
reaction,  releasing  additional  energy  as  A1  oxides  are  formed. 

The  cracking  and  leaking  of  A1  onto  the  surface  of  the  Ni  shell  (figure  14)  is  very  different  from 
what  Delogu  (18)  observed  using  a  semi-empirical  tight-binding  (TB)  potential.  Another 
difference  from  the  work  by  Delogu  is  the  choice  of  A1  core  radius.  For  both  of  the  shell  cases 
here,  the  core  has  a  radius  of  2.8  nm,  this  is  in  contrast  to  the  core  radius  of  3.0  and  2.0  nm  for 
the  1.0-  and  2.0-nm  shell  cases,  respectively,  used  by  Delogu  (18).  Since  the  1-nm  shell  case  is 
the  most  interesting  and  the  radii  are  very  close  (2.8  vs.  3.0  nm),  the  choice  of  core  size  is  not 
expected  to  have  had  an  appreciable  effect  on  the  results  for  this  case. 


Figure  14.  Cross  sections  of  the  Ni-coated  A1  nanopartiele  simulation  model  just  after  melting  of  the  A1 
core  for  the  (a)  1-nm-thick  and  (b)  2-nm-thick  Ni  shells. 


7,3  Analysis  of  Alloy  Formation 

During  the  sintering  of  Ni  and  A1  nanoparticles,  numerous  local  atomic  ratios  may  be  identified 
from  the  Ni-Al  binary  system  (32-34).  The  data  presented  here  are  computed  from  the 
simulation  results  by  calculating  the  local  atomic  ratios.  These  ratios  are  computed  in  order  to 
analyze  the  spatial  and  temporal  distribution  of  alloys  in  the  simulation.  The  most  commonly 
observed  alloys  for  the  Ni-Al  binary  system  include  AlgNia,  AfiNi,  Al3Ni2,  NiAl,  NijAfi,  and 
NisAl  (55)  corresponding  to  Ni  fractions  of  0.18,  0.25,  0.40,  0.50,  0.625,  and  0.75,  respectively. 
A  comprehensive  analysis  of  the  measured  atomic  ratios  in  solution  during  sintering  of  separate 
nanoparticles  is  summarized  here.  Figures  15a  and  15b  show  the  normalized  amount  of 
computed  solution  ratios  by  Ni  fraction  through  stage  2  and  during  stage  3,  respectively. 
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Figure  15.  The  results  here  show  the  fraction  of  Ni  in  each  alloy  at  stages  (a)  1  and  (b)  2  versus  fraction  of  all  atoms 
in  each  alloy,  i.e.,  a  data  point  at  (0.25,  0.1)  corresponds  to  10%  of  all  atoms  in  the  system  belong  to 
AljNi. 


During  stage  1,  from  0  to  50  ps  for  sintering  of  10-nm-diameter  partieles,  most  of  the  atoms 
remain  in  the  pure  metal  phase,  with  only  a  small  fraction  of  atoms,  less  than  5%,  forming  Ni 
rich  solutions  such  as  NiaAl.  In  figure  6a,  at  50  ps,  an  observed  peak  at  0.75  fraction  Ni  (NisAl) 
is  likely  due  to  the  shorter  bond  length  of  Ni-Ni,  such  that  near  the  contact  interface  each  A1 
atom  will  bond  with  multiple  Ni  atoms.  During  stage  2,  the  formation  of  A1  rich  liquid  solutions 
occurs  at  a  rapid  rate,  quickly  depleting  the  atomic  fraction  of  pure  A1  in  the  system  by  the  end  of 
stage  2  (500  ps).  In  figure  1 6a,  at  500  ps,  a  small  peak  between  0.2  and  0.3  fraction  Ni  is 
observed,  corresponding  to  A^Ni.  At  the  end  of  stage  2,  the  largest  portion  of  any  material  is  still 
pure  Ni  (figure  16b).  The  initial  formation  of  A1  rich  solutions  is  supported  by  experimental 
observations  (32,  34).  The  reaction  process  that  produces  these  alloys  is  likely  explained  by  Ni 
atoms  diffusing  into  the  liquid  A1  region  forming  A1  rich  clusters  of  AfiNi.  These  clusters  then 
diffuse  through  the  A1  melt,  eventually  forming  the  eutectic  mixture.  Figure  16b  supports  this 
process,  because  small  clusters  of  A1  rich  alloys  are  clearly  evident  in  the  liquid  Al,  with  some 
initial  coalescence  of  clusters  to  obtain  the  eutectic  atomic  ratio. 


Figure  16.  Plot  showing  spatial  distribution  of  Ni  fractions  during  stages  (a)  1  and  (b)  2  of  the  sintering  process. 
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The  lines  in  figure  17b  for  stage  3  indicate  a  rapid  reduction  in  the  amount  of  pure  Ni  and  the 
corresponding  formation  of  NiAl.  During  stages  1  and  2,  as  observed  in  figure  3,  the  temperature 
of  the  system  is  rising  monotonically.  Near  the  start  of  stage  3,  the  system  temperature  has 
reached  the  melting  point  of  Ni.  At  700  ps,  shortly  after  the  Ni  nanoparticle  has  completely 
melted,  two  peaks  are  present  in  figure  17a  near  40%  Ni  and  50%  Ni,  corresponding  to  A^Nii 
and  NiAl.  In  figure  17a,  shortly  after  melting  of  the  Ni  nanoparticle,  there  is  a  Ni  rich  region 
visible.  This  region  quickly  disappears  as  the  SHS  reaction  proceeds  rapidly  to  form  the  eutectic 
alloy.  After  nearly  complete  sintering  at  950  ps,  the  majority  of  atoms  are  in  a  NiAl  liquid 
solution  with  some  temporal  formation  of  other  alloys.  This  brief  analysis  has  demonstrated  that 
the  various  solution  ratios  are  an  interesting  problem  that  warrants  closer  thermodynamic 
analysis. 


Figure  17.  Plot  showing  spatial  distribution  of  Ni  fractions  at  the  (a)  beginning  and  (b)  end  of  stage  3  of  the 
sintering  process. 


8.  Conclusions 


We  have  analyzed  two  model  systems  for  the  energetic  reaction  of  Ni  and  Al.  In  the  first  case, 
we  considered  the  coalescing  and  sintering  of  separate  nanoparticles  and  found  that  the  energy 
release  from  the  change  in  surface  area  was  only  significant  for  small,  less  that  10-nm-diameter, 
nanoparticles.  These  separated  nanoparticle  reaction  simulations  and  thermodynamic  analyses 
showed  that  the  reaction  time  decreased  and  the  adiabatic  reaction  temperature  increased  with 
decreasing  nanoparticle  sizes.  This  may  be  important  for  applications  where  high  energy  release 
rates  are  desired.  The  simulation  data  closely  matched  a  classical  thermodynamic  analysis. 

In  the  second  part  of  this  work,  we  considered  the  sintering  of  Al-coated  Ni  nanoparticles  and 
Ni-coated  Al  nanoparticles  as  a  model  material  system  for  nanoparticles  embedded  in  a  matrix  of 
the  other  metal.  This  work  revealed  that  the  reaction  time  was  again  inversely  related  to 
nanoparticle  size,  but  the  adiabatic  temperature  decreased  with  decreasing  nanoparticle  size. 
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Mechanically,  the  Al-coated  Ni  nanoparticle  system  was  a  model  system  for  a  lightweight  A1 
matrix  with  embedded  Ni  nanoparticles,  a  system  with  relatively  high  strength  compared  to  a 
loosely  bonded  powder  of  A1  and  Ni  nanoparticles.  This  A1  matrix  system  could  be  used  in 
systems  where  mechanical  strength  is  important  in  addition  to  energy  release  from  kinetic 
sintering  of  the  Ni  and  A1  atoms.  In  the  Ni-coated  A1  nanoparticle  system,  we  investigated 
possible  rupture  and  fragmentation  of  the  Ni  shell  but  were  unable  to  observe  any  fragmentation. 
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List  of  Symbols,  Abbreviations  and  Acronyms 


A1 

aluminum 

ARO 

U.S.  Army  Research  Office 

EAM 

embedded  atom  method 

LAMMPS 

Large-scale  Atomic/Molecular  Massively  Parallel  Simulator 

MD 

molecular  dynamics 

MSRC 

U.S.  Army  Major  Shared  Resource  Center 

Ni 

nickel 

NiAl 

nickel  aluminide 

NIST 

National  Institute  for  Standards  Technology 

SHS 

self-propagating  high-temperature  synthesis 

TB 

tight-binding 
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